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Searches for lepton flavor violating decays £^7, ^ ^^P*^ (where 

£~ = e^,fi~, and = 7T^,r],r]') at B-Factories: Status and Combinations 

Swagato Banerjee^* 

^Department of Physics and Astronomy, University of Victoria, 
P.O. Box 3055, Victoria, British Columbia, V8W 3P6 Canada. 

The present experimental status of searches for lepton flavor violating decays of the r lepton to a lighter mass 
lepton and a photon or a pseudoscalar meson is presented. The results obtained are from e^e~ annihilation data 
collected at a center-of-mass energy near 10.58 GeV by the BABAR and the Belle detectors. In order to benefit 
from larger available datasets, the first BABAR— Belle combinations of upper limits on the branching fractions of 
lepton flavor violating r-decays are obtained. 



1. Introduction 

The recent discovery of neutrino mixing [T] im- 
plies that lepton flavor violation (LFV) occurs. 
LFV in charged decays have not yet been ob- 
served, although they have long been identified 
as an unambiguous signature of new physics. 

The Standard Model (SM) extended to include 
finite mass differences from neutrino mixing im- 
ply T — /i mixing is generated at the one-loop level 
and thus suppressed by a factor of (m^/m^)^, 
which is many orders of magnitude below the 
experimental sensitivity. However, many new 
theories, for example SM with additional heavy 
right-handed Majorana neutrinos or with left- 
handed and right-handed neutral isosinglets [3]; 
mSUGRA models with right handed neutrinos 
introduced via the seesaw mechanism [4|5j ; su- 
persymmetric models with Higgs exchange |6|7) 
or SO(IO) symmetry |8|9j: or technicolor mod- 
els with non- universal Z' exchange fW], allow for 
decays, where £~ = e^,/i^, at the 
level of ~ 0(10-1° - 10"^). 

Large neutrino mixing [1] could potentially in- 
duce large mixing between the supersymmetric 
partners of the leptons. While some scenario's 
predict higher rates for — > /^*7 decays, oth- 
ers, for example with inverted hierarchy rather 



than normal hierarchy of masses of supersymmet- 
ric partners of leptons, predict [5] higher rates for 
e^7 decays. 

Semi-leptonic neutrinoless decays involving 
pseudoscalar mesons like 



where 



*To appear in the proceedings of g*** International Work- 
shop on Tau-Lepton Physics (Tau06) held at Pisa, Italy, 
from 19"* to 22"'' September 2006. 



P° = 7r°, 77,77', are likely candidates for LFV in 
supersymmetric models [llj . for example, aris- 
ing out of exchange of CP-odd pseudoscalar neu- 
tral Higgs boson, which are further enhanced by 
color factors associated with these decays. The 
large coupling of Higgs at ss vertex enhances fi- 
nal state containing the 77 meson, giving a pre- 
diction of B{t^ fJ-^v) ■ 13{t^ — + /i''^7) = 
8.4 : 1.5 jT2]. Some models with heavy Dirac 
neutrinos |13|14j . two Higgs doublet models, R- 
parity violating supersymmetric models, and fla- 
vor changing Z' models with non-universal cou- 
plings [H] allow for observable parameter space 
of new physics [16j , while respecting the existing 
experimental bounds at the level of ~ 0(10"^). 

Present generation e+e" B-factories operat- 
ing around the T(4S) resonance also serve as r- 
factories, because the production cross-sections of 
'^bb ■ '^r+T~ = 1-1 • 0.9 nb are comparable to 
each other at center-of-mass (CM) energy near 
10.58 GeV. Recent results from the searches of 
i^j/P'^ are presented here, along with a 
first attempt at combining the upper limits from 
the BaBAR and the Belle experiments. 
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2. B-Factories: Status 

In e~^e~ annihilations, £^^/P^ decays 

have 2 characteristic features: the energy of tau- 
daughters is close to the half the CM energy, and 
the total invariant mass of the daughters is close 
to the mass of r-lepton. Sometimes, the mass 
of the daughters are calculated with total energy 
constrained to be half the CM energy, in order to 
improve upon the resolution of the mass variable. 

The general strategy for searches of such de- 
cays is to define a signal region in the energy- 
mass plane of the r-daughters (typically 2a re- 
gion), and optimize upon a set of selection criteria 
to reduce the contribution to background events 
from well-known SM processes. The principal 
sources of background are radiative di-muon or 
Bhabha events, along with some irreducible con- 
tribution from r-pair events having hard ISR with 

— > i^iyp or TT^TT^v decays for 

or £^P^ searches, respectively. 

After applying the selection criteria, the sig- 
nal efficiency (e) and expected background events 
(TVbkg) are estimated inside the signal region. Fi- 
nally, the number of observed events (A'obs) are 
counted (or fitted) from the data. As yet, no 
signal has been observed in any of the searches. 
So, an upper limit at 90% confidence level (C.L.) 
is obtained as: B'^^ = N^JiNrS), where TV^J^is 
the 90% confidence level (C.L.) upper limit on the 
number of signal events inside the signal box, and 
Nt = 2C(7^+T— is the number of r-decays studied 
from a sample of data with a luminosity of £. 

For the £^P^ searches, the pseudoscalar 

mesons are reconstructed in the following decay 
modes: tt'^ 77 for —> i^n^, 77 — > 77 
and r] tt+tt^tt^ {n^ — > 77) for £^V: 
r( TT^n^r] {rj —> 77) and 7/ — s- p'^7 for 
— > £^?]' . In these searches, the signal efficiency 
(e) in the above formula for is the product 
of the branching fraction (B) and reconstruction 
efficiency (^) of the P" decay mode under consid- 
eration. To obtain a combined upper limit with 
?7 and 77' decays, the observed and expected back- 
ground events and the signal efficiencies are added 
using BS, = {Bi X ^1 -I- X ^2), where Bi, B2 are 
the respective branching fractions and ^1, ^2 are 
the corresponding reconstruction efficiencies. 



Searches for LFV in t± ^ ^±7 and t± £^P° 
decays have been performed with C — 232.2 fb^^ 
and 339 fb"^ by the BaBAR experiment |2ll7|18j . 
and with C = 535 fb"^ and 401 fb"^ by the Belle 
experiment ^19.20] . respectively. The results for 
these searches are summarized in Table 1, where 
B{tt° 77) = (98.80 ± 0.03)%, B{ri 77) = 
(39.38 ± 0.26)%, 6(77 -> 7r+7r-7r°,7rO ^ 77) = 
(22.43 ± 0.40)%, B{ri' n+n-rj,T] 77) = 
(17.52 ± 0.56)%, and 6(77' ^ p°j) = (29.40 ± 
0.90)% have been used [21 for combining modes. 

3. B-Factories: Combinations 

The results from the BaBAR^ and the Belle ex- 
periments are combined by adding the number of 
observed and expected background events along 
with their symmetrized errors, respectively, in- 
side the signal region. The signal efficiencies are 
averaged using weights than account for the inte- 
grated luminosity for the corresponding searches. 
This is the most conservative approach given that 
the signal and background PDFs are different for 
similar searches by different experiments, because 
of differing resolutions on reconstructed mass and 
energy among the experiments. Thus, in this pro- 
cedure, the distribution of both the signal and 
the background events inside the signal region are 
taken to be uniform, and only the rates of signal 
and backgrounds events contribute to the combi- 
nation. 

Combined frequentist limits are calculated in- 
cluding all uncertainties using the technique of 
Cousins and Highland ^22 following the imple- 
mentation of Barlow [23]. In this technique, 10^ 
MC samples are generated according to a Pois- 
son distribution with mean (s -I- 6), where the 
background (6) and signal (s) are each drawn 
randomly from Gaussian distributions describ- 
ing their respective PDFs. For the signal Gaus- 
sian, the mean and the standard deviation are 
2CaTTB\jh£ and its uncertainty, with the latter 
determined by propagation of errors. The branch- 
ing ratio (.Bul) is varied until we find a value 
for which 10% of the samples yield a number of 
events less than that observed in the combined 
data. At 90% C.L. this procedure gives the com- 
bined upper limit observed for each decay mode. 
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Table 1 

The luminosities, the efBciencies, the number of expected and observed background events inside the 
signal region, and the expected and observed upper limits on the branching fractions at 90% C.L.{B^j^) 
for the searches from the BaBAR |2I17I18| and the Belle [19l20j experiments and their combinations as 
described in the text. An empty cell indicates that the respective value has not been published. 
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In order to calculate the expected upper limit in 
absence of a signal, the upper limit for each guess 
of the observed events is weighted by a Poisson 
distribution centered around the number of ex- 
pected background events and averaged. 

The results of these combinations for the ex- 
pected and the observed upper limits at 90% C.L. 
(Sul) ^re summarized in Table 1. In the case 
of — > /i*7, the observed limit is significantly 
lower than the expected limit, because both ex- 
periments observe less events than expected [24] , 

4. B-Factories: Projections 

Depending upon the nature of backgrounds 
contributing to a given search, two scenarios can 
be envisioned while extrapolating to higher lu- 
minosities. In the first scenario, there is no sig- 
nificant background contribution, and the selec- 
tion procedure may be optimized to give 0(1) 
expected background events, and so, Bfji^ oc l/£. 
In the other scenario, there is an irreducible back- 
ground contribution, so that A^^^ ^/Ndbs grows 
as VZ, and so, cx I/VZ- 

Taking as baseline the expected upper limit of 
B{t^ M^7) = 1-2 X 10~^ from search by the 
BaBAR experiment with C = 232.2 fb^^, one can 
expect upper limits of S(r^ M^t) — 1-4 x 10"^ 
or 4.1 X 10~^ in scenarios with background free 
or background dominated search, respectively, in 
a couple of years when both the BABAR and the 
Belle experiments have accumulated total data 
with luminosity of £ = 2ab^^. If no signal is 
found, the observed limit may be worse or better 
than the expected limit, depending upon if the 
observed number of data events is more or less 
than the expected number of background events, 
as evident from the trend in Table 1. 

5. Constraints on new physics 

Improved limits on LFV in r decays give several 
interesting constraints on new physics. Two such 
constraints are discussed here: one from B{t^ — > 
/i^7) and one from B{t^ — > iJ.^ri). 

In supersymmetric models with SU(5) symme- 
try for Grand Unification Theories, flavor chang- 
ing right-handed currents induce correlations be- 



tween possibly enhanced rates for r — > 7 de- 
cays and CP-asymmetry in b-s penguins |25j . The 
HFAG summer 2006 average [25 of BaBAR [27] 
and Belle ^28j measurements of S{B <j)K^) = 
(0.39 ± 0.18) is already very promising. How- 
ever, an observed upper limit of B{t^ M^7) 
< 1.6 xl0~* from this combination of results 
from the BaBAR and the Belle experiments with 
a combined C — 767.2 fb^^ is more constrain- 
ing, as shown in Figure 1. The contours are 
shown for fixed values of universal gaugino mass 
rUg = 400, 600, 800, and 1000 GeV, with uni- 
versal scalar mass toq S [200, 1000] GeV, the ra- 
tio of the vacuum expectation values of the two 
Higgs doublets tan/? = 10, the trilinear coupling 
Aq = 0, the mass of the right-handed neutrino 
"^i/jj = 5 X 10^"* GeV, and the mass of left-handed 
T-neutrino m^^ = 5 x 10~^eV (for description of 
other parameters see Figure 2 in Reference [25]). 

Mixing between left-handed smuons and staus 
allows one to translate the upper limit on B{t^ 
H^rj) to an exclusion plot in the tan/3 vs. 
plane, where niA is the mass of the CP-odd pseu- 
doscalar neutral Higgs boson. The excluded re- 
gions at 95% C.L. from this combination of re- 
sults from the BABAR and the Belle experiments 
for — > /i^77 search (< 6.8 x 10^*) with right- 
handed neutrino mass = 10^** GeV introduced via 
the seesaw mechanism in supersymmetric mod- 
els [12] are shown in Figure 1. Also are shown 
the regions excluded by the BABAR T8j and the 
Belle [20] experiments individually. These re- 
sults are competitive with those obtained from 
the direct searches for Higgs bb, t^t~ de- 
cays by the CDF [29] and the DO [30] experi- 
ments, and complementary to the region excluded 
by the LEP experiments with a top quark mass 
of 174.3 GeV [31], for two common scenarios of 
stop-mixing benchmark models [32]: m™'^'^ and 
no-mixing obtained with the Higgs mass param- 
eter ^ = —200 GeV shown by darker and lighter 
shaded regions respectively. 

6. Summary 

An improvement of five order of magnitude in 
the upper limits for B{t^ 6^7) and B{t^ 
/i*7) over the last twenty-five years [21j is shown 
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in Figures 3 and 4, respectively. The next five 
years promises to be the most interesting phase in 
this evohition, when experiments approach closer 
to the predictions of different theoretical models. 
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Figure 3. Evolution of experimental bounds of 
B{t^ 6^7) and some theoretical predictions. 



Figure 1. S(r± ^i^-f) vs. S{B (j)K°) con- 
tours as function of jjig, mo (see text). 
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Figure 4. Evolution of experimental bounds of 
B{t^ H^j) and some theoretical predictions. 



Figure 2. Excluded regions in tan/3 vs. plane 
from search for t"^ fx^rj decays (see text). 
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